arXiv:1509.04341v2 [quant-ph] 25 May 2016 


Deriving Entropic Inequalities for Two Coupled Superconducting Circuits 


Evgenii Glushkov^’^*, Anastasiia Glushkova^ and V. I. Man’ko^’^ 


^Moscow Institute of Physics and Technology 
Institutskii per. 9, Dolgoprudnii, Moscow Region If 1700, Russia 
‘^National University of Science and Technology ’’MISiS” 

Leninskii pr., f, Moscow 119991, Russia 
^Lebedev Physical Institute, Russian Academy of Sciences 
Leninskii pr., 53, Moscow 119991, Russia 
* Corresponding author e-mail: eugene.glushkov@gmail.com 

Abstract 

We discuss the known construction of two interacting superconducting circuits, based on Josephson 
junctions, that can be precisely engineered and easily controlled. In particular, we use the parametric 
excitation of two circuits, realized by an instant change of the qubit coupling, to study entropic 
and information properties of the density matrix of the composite system. The density matrix is 
obtained from the initial thermal state and is then analyzed in the approximation of small perturbation 
parameter and sufficiently low temperature. We also check the subadditivity condition for this system 
both for von Neumann and deformed entropies and look at the dependance of mutual information on 
the temperature of the system. Finally, we discuss the applicability of this approach to describe such 
system of two coupled superconducting qubits as harmonic oscillators with limited Hilbert space. 
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1 Introduction 


The idea to use the Josephson junction to engineer superconducting circuits, where current and voltage 
are considered analogs of position and momentum in a parametric oscillator , is now employed to study 
the properties of qubits, associated with the state of such circuits. The non-stationary quantum states of 
current and voltage in such circuits have been extensively studied [2j|^ for more than two decades and 
enormous progress in experimental realization of these systems has been made, which gave rise to 
the whole new area of research called quantum information processing. The recent works [8-11 study 
the information properties of qubits as the resource for future quantum technologies. The particular 


properties of the multi-qubit states to be analyzed are entropic and information inequalities 13 
which serve as the basis for quantum information processing. 

In this work we aim at studying the subadditivity condition for von Neumann and q-entropies of the 


bipartite quantum system 14-19 , on the specific example of two coupled superconducting circuits 20 


To study entropic and information properties of the composite system we use the parametric excitation 
of two circuits, realized by an instant change of the qubit coupling. Initially the system is considered 
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to be in a thermal state with a corresponding thermal density matrix, which is then analyzed in the 
approximation of small perturbation parameter and sufficiently low temperature. We also look at the 
dependance of mutual information on temperature of the system and discuss the applicability of our 
approach to describe two coupled superconducting qubits as harmonic oscillators with limited Hilbert 
space. 


2 Theoretical model 

We start with a simple model of two interacting harmonic oscillators in a thermal state. As we would 
like to consider only the first two levels in each of the circuits, we are working in the limit T —)• 0. This 
approximation and its applicability will be more thoroughly discussed in Section 



Figure 1: Schematic representation of two superconducting circuits, modelled as harmonic oscillators, 
coupled by a mutual inductance. 

The Hamiltonian of the system of two coupled resonant circuits, depicted in the Fig. [^can be written 
in the following simple form: 


H = H 1 + H 2 + V, 


( 1 ) 


where first two terms correspond to two independent LC-circuits and the third term dehnes coupling 
between them: 


= ^ + H = Li2/i/ 2, i = l,2 

For the convenience and due to the duality between mechanical oscillators and LC-circuits, we now 
introduce canonical position and momentum operators, which we will be using hereinafter, through the 
following change of variables [21] : 

% = Pi = = ifiSjk, j, k = l,2 
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Using this substitution we can rewrite Eq. in the following form: 


H = 


Pi 

2m 


+ 


mujrx 


, rnulxl 
^ 2m^ 2 


+ gmxiX2UJiUJ2 


( 2 ) 


where g = Li 2 {LiL 2 )~^/‘^ is the qubit coupling constant and are the eigenfrequencies of the circuits. 
For simplicity, here and later we assume that m = 1, wi = 1, W2 = A. 

Next we would like to get rid of the cross term and diagonalize the Hamiltonian [22] . To do that 
we apply the following rotation by an angle (f), using the theorem that a quadratic form is reduced to a 
diagonal form by an orthogonal transformation: 



After this transformation, we obtain the following Hamiltonian: 


p\ p2 cos^ (l)x'i + sin^ 4>X2 


/2 


H=^+^+ 
2 2 


+ cos (j) sin (l)X 2 x'i+ 


+ 


A^(sin^ 4>x[ + cos^ 0x2) 


— A^ sin (j) cos (f)X 2 Xi+ 


+ gX{— cos 4> sin (j)Xi + cos </> sin 4>X2) + gX{cos^ (p ~ sm (p)x 2 Xi (4) 


Equating to zero the cross terms from Eq. we obtain the system of two independent harmonic oscilla¬ 
tors: 


p'l Vt\x'i p '2 Xlox' 

H = — -\ -+ — + 

2 2 2 2 


2 

2-^ 2 


(5) 


where the new resonant frequencies and the rotation angle in the approximation of small angles are: 

« 1 — 2gXc/) -I- A^^^, 142 ~ + A^ + 25'A(/), p ^ gX/ (A^ — 1) . (6) 


3 Density Matrix 


As our oscillators are in the thermal state in the low temperature limit the density matrix of the system, 
reduced to the 4x4 subspace, will take the following form: 


/ exp 


^no 

T 


1 

W) 


\ 


0 

0 

0 


0 0 

-thi 

exp T 0 

„ -hm 

0 exp T 

0 0 


0 ^ 

0 

0 

_£ii , 
exp T j 


(7) 


where Enm = f4i(n + 1/2) + r42(m + 1/2) and .Z’(T) = Zi{T) ■ Z 2 {T) = (4sinh ^ sinh ^) ^. 

However, if we want to return to our initial system with two coupled resonant circuits, we have to 
perform transformation, which decomposes the old eigenstates in the new rotated basis: 


p = U-^p U 


( 8 ) 
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In the basis of the harmonic oscillators eigenstates we denote the coefficients of the transformation 
Unmn'm' and Calculate them using eigenfunctions of the harmonic oscillators: 


Ur, 


x-^ X2 X2 

VT 7-00 \/fLiL 2 V 2 ’" 2 ™n!m! \/2^'2'^'n'lm'l ^ I ^ 




(9) 


where: 


I = 




(jj — A, Li — 




and Hi are the corresponding Hermite polynomials. 

After calculating the matrix elements of the transformation Unmn'm' (see Appendix for the details) 
we obtain the hnal density matrix to work with from the Eq. which can be seen in the Fig. at 
various temperatures of the system. One can see, that the approximation used is valid only for very 
low temperatures, below 100 mK, which is standard for superconducting qubits. At higher temperatures 
the populations of the higher energy levels become non-negligible, thus demanding to take into account 
larger subspace, than the 4x4 matrix we are discussing here. 



Figure 2: (color online) Density matrices, calculated from Eq. Q, at different temperatures. 
Off-diagonal elements are nonzero, but not visible due to the scale difference. 


4 Calculating entropy 

Using the obtained density matrix, we can calculate the entropical properties of the system. The main 
properties we will be looking at are entropies of the bipartite system and its subsystems (single qubits) 
and the mutual information, which can be deduced from the subadditivity condition. 

The two types of entropies we are discussing in this work are the von Neumann entropy [S = —Trp In p) 
and the Tsallis entropy, which equals the von Neumann entropy in the limit of (7 —)• 1: 

Sj =-Trplnqp, (10) 

where the q-logarithm is defined as follows: 
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ln9>o P 


[lnp,g = 1 


As the matrix elements of the density matrix depend on temperature of the system, we calculated the 
Eq. [^for various temperatures and displayed the results of these calculations in Fig. The temperature 
here and later is chosen to be in units of frequency, relative to the eigenfrequency of the qubits oj. 



0 0.05 0.1 0.15 0.2 

Temperature (in dimensionless units) 


Figure 3: (color online) Dependance of the q-entropies on temperature calculated from Eq. 


10 


One can see, that the q-entropies in Fig. [^are grouped in such a way, that the entropies with q <1 
lie higher than the von Neumann entropy {q = 1), while the entropies with g > 1 lie lower. The meaning 
of this behavior can be understood in the following way. The higher is the value of q, the higher is the 
influence of the biggest terms in the distribution on entropy and the more deterministic is the behavior 
of the system. As the main elements of our density matrix are the four diagonal ones, higher q-values 
lead to lower value of q-entropy for our matrix. 
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5 Verifying subadditivity condition 

To check the subadditivity condition for the bipartite system, consisting of two resonant circuits, we 
divide the calculated density matrix from Eq. [^into density matrices of the subsystems as follows: 

Pll + P 22 Pl 3 + / 524 \ _ ( Pn + P33 P 12 + P34 

’ ^2 = 

P31 + P42 P33 + Pii J \P21 + PA3 P22 + PAA 

And the subadditivity condition for the corresponding q-entropies reads: 

S^{p)<S^{pi) + S^{p2) ( 12 ) 

One can also deduce the value of mutual information out of it: 

I = Sjip) - (pi) - (P2) > 0 (13) 





Figure 4: (color online) Mutual information related to temperature. 
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Using Eq. 13 we calculate the value of mutual information for various temperatures of the system 
and present them in Fig. The fact that the value of mutual information stays positive in the whole 
temperature range shows that the subadditivity condition holds for our system. 


6 Approximation applicability 

In this section we discuss the applicability of the approach we used to approximate the system of two 
harmonic oscillators by two qubits in the limit of low temperature. We look at the density matrix of our 
system and calculate the purity parameter /r = Tr , first for the reduced 4x4 matrix (I) and then for 
the remaining part of the density matrix (II). The result of these calculations is shown in the Fig. 5a 
in relation to temperature of the system. We also plot the sum of non-diagonal elements of the reduced 
4x4 matrix to see the region, where it’s bigger that the purity error. 


Additionally in the Fig. 5b we plot the purity of the 4x4 density matrix versus temperature to show 
the applicability of our approach. One can see that the 4x4 subspace approximation is valid below 
T PS 0.2, where p ps 1 . 

Based on these calculations, we believe that the approach used is valid for temperatures up to 100 mK, 
or 0.2 in relative units. So, for all the calculations with the reduced 4x4 density matrix, we normalize 
its elements to get p = 1. 



(a) Relation of purities. 

Figure 5: (color online) Applicability limits, (a) T 
approximation to the purity of the remaining part 
of the 4-level approximation on temperature. 



(b) Purity for different temperatures. 

inverse relation of purity of the 4-level 

the density matrix, (b) The dependance of purity 


7 Summary and conclusions 

We have analyzed the system of two coupled superconducting circuits, modeled as two interacting har¬ 
monic oscillators in the low-temperature limit, and derived its density matrix in the approximation of 
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small perturbation. We have further checked that the calculated density matrix of the bipartite sys¬ 
tem satisfies the entropic inequalities both for the von Neumann and Tsailis entropy and looked at the 
dependance of mutual information on the temperature of the system. Finally, we evaluated the purity 
parameter of the system and verified the applicability of our approach to describe the system of two 
coupled superconducting qubits as harmonic oscillators with limited Hilbert space. 

As the next step we would like to compare our findings with experimental data and generalize our 
approach for arbitrary rotation angle cj) and for larger number of qubits, which we will do in future 
publications. 
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Appendix 

Here we describe in more detail how the matrix elements were calculated. In the basis of the eigenstates 
|m, n) we denote the coefficients of decomposition Unmn'm' and calculate the matrix elements of the U in 
the eigenvalues basis of the Hamiltonian (Eock basis). 

Eor the transition from an old basis to a new one \ek) —?• |e~fc) one can write the following decomposition: 


[Cfc) — ^ ^ \Cm) ^ ^ {(^n\ Ukm \Sm) — {^n\ 1 \^k) — 
m m 

— ^ ^ Ukm (^nl i [Cjn) — ^ ^ U^mdum — U^n — 


k{x)dx (14) 


So, for the density matrix we obtain: 


Pnmn'm' — {rim\ p\n m) 


1 

W) 


inm e 


-pH 


! I\ 

n m ) 


(15) 


Using the harmonic oscillator eigenfunctions: '^n{x) 




Hn{x), where oj 


1 , we get: 
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u„. 


_4_4 

r°° g 2 g 2 2 


dxi 


dX2 


TT 


y/l \/LiL 2 V2'^2'^n\m\ \/2^'2^'n'\m'\ 

X2' 


Hn (Xl) Hm ( y) i^n' ) = 


1 

TT 


O O ^'2 ^/2 

^2 -r.2 ^2 


1 


1 


-oo \/ L 1 L 2 \J2^2 ™n!m! 


Hn (Xl) Hm{f]Hn 


X2 


COS 0X1 + sin 0X2 




— sin 0xi + cos 0 X 2 
L 2 


roo g-^-5,4 


TT 


_ gn(xi) Hm{^) 

VI V 2 ^' 2 ™'n'!m'! V 2 ” 2 ™n!m! 




cos 0xi + sin 0 X 2 




— sin 0xi + cos 0 X 2 
Lo 


(cos03;X+Bin<j>rg2) ^^2 gjj^2 0_|_(.os^ tf>—2g\ s'm (f) COS (j)) — - — '^^2) (^gjp2 pQg2 0_|_2(j,A s'm tf> COS tf>) 


VLIL 2 


-dxidx2 


where I = ^ ^\,uj = \Li = ^H^L 2 = ^ 

sin^ 0 + cos^ 0 — 2g'A sin 0 cos 0 
^^2 = sin^ 0 + A^ cos^ 0 + 2g\ sin 0 cos 0 


In the approximation of small angles (cos 0 ~ 1 and sin 0 ps 0) last terms are: 


Hn 


COS 0xi + sin 0X2 
. ~Li 

n\ ^ + 1 - 2c/A0 

~ 0^ + A^ + 2gX(f) 


Hn 


— sin 0xi + cos 0 X 2 
L 2 


Hn 


Xl + 0X2 


Hn 


X2 - 0Xl 
L 2 


Here we denote K = VHW ~ (A^(A^0^ + 1 — 25 rA 0 )( 0 ^ + A^ + 2grA0)) and use Ho{x) = 1, Hi{x) 

2x. With these simplifications we obtain the following expression for the matrix elements: 


Un 


— f Q-Vi^l+4'‘^^2+l)xl+{<j>^li-<t>^2)xiX2 + ^{<j>^Ul+^l2+>y)x2) (^Xl) Hn 

ttK J_^ 

dxi dx2 


0X1+0X2 0 „ 0X2 -0X1 

• iirt/ I - I ^ra' I - / _ _ 

Li J \ L 2 J V2'^2^n\m\ V2«'2™'n'!m'! 

Using this equation, we now calculate the matrix elements: 


t^oooo — 

Ktt 


roo 2,2 ^2 q^^'2 


g 2 ^ 


2 dxidx2 


(16) 


(17) 
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V2{cj)^ + X^ + 2gX<P)^/^ ^ ^ ^ 

Uoooi = -TTT- / (®2 - (l)Xi)e 2 W 2 2 dxidx 2 


2 2 f2 /2 

x-^ X2 \I 2 X 2 


ttK 


( 18 ) 


[/ooio =-- / (a^i + (t>X 2 )e 2^2 2 dxidx 2 

VTA ./_^ 


(19) 


t^ooil — 


^ a ^2 ri 2 a; 2 ^ 


TT 


^/AK3 y_. 


(xi + (^X2)(a:2 — i;/>xi)e 2^2 2 dxidx2 


( 20 ) 


t^oioo — 


Ktt 


X2e 2 ^ 


2 2 i2 i2 

x^ !.i\x-^ 1223^2 


~ dxidx2 


( 21 ) 


t^oioi — 


2 2 o 22 o /2 

X-^ X 2 ^l-^X-^ \I 2 X 2 


TTK^{X^(fP‘ + 1 — 2gX(j)Y/^ y_c 


X 2 {x 2 — xicj))e 2 w 2 2 dxidx 2 


( 22 ) 


2 2 o /2 /2 

3^2 s2-]^a:2 1223:2 


= nKH^^ + A 2 + 2,A^)i/4 U ^ 


(23) 


t^oiii — 


2^/2 

ttK^ 


2 2 o ^'2 o ^'2 


a:j^ a :2 ii 2^'2 


0-2^ 


X2{X2 — (t)Xi){xi -\-(l)X2)e 2^2 2 dxidx2 


(24) 


t^iooo = -r^ / Xie“^“2i2 2 ^dxidx 2 

Ktt ,/_^ 


(25) 


2 2 /2 o /2 

^2 i2]^x^ 1223:2 


2((/.2 + A2 + 2(7A(/.)V4 /•oo ^ 

t/1001 =- ^7 - / Xl(x2-0Xl)e 2 W 2 2 dxidx2 

Ktt ,/_^ 


(26) 


t^ioio — 


2 (A 202 + 1 - 2 c/A(/>)^/^ 


Ktt 


/ OO 

Xl{xi + (j)X2) 

-OO 


2 2 r-, /2 n, /2 

^2 12-j^x^ 1223:2 


e 2 ^ 2 2 dxidx 2 


(27) 


t^ioii — 


2\/2 r 
■VXK^ J-c 


2 2 /2 o /2 

X_ ^2 ^1^1 ^ 23:2 


xi{x2 — (pxi){xi + (pX2)e 2^2 2 dxidx2 


(28) 
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t^lioo — 


2 ^/A 

Ktt 


2 2 /2 /2 
X 2 \I 2 X 2 

xiX2e~^~^ ^ ^ dxidx2 


(29) 


t^iioi — 


2 V 2 


TTK^{X‘^(j)‘^ + 1 — 2gX(l)Y/^ J_c 


XiX 2 {x 2 - <t>xi)e 2 2 ! 


2 2 o /2 o /2 

1 X2 CLix-^ 02^2 

2 2 dxidx2 


(30) 


t^iiio = 


2\/2 


/ oo xf x| njxf 

xiX2(xi + 4>X2)e~^ ^ 2 2 dxidx 2 

-00 


7rK^{(p‘^ + + 25 A(/))^/‘^ j_ 


(31) 


t /1111 = 




f 


2 2 o /2 o /2 

2 2^2 SZ 2 a ^2 


XiX2(x2 — (/>Xi)(xi + (/>X2)e 2^2 2 dxidx2 


(32) 


Calculation of matrix elements 

To calculate the integrals we use the following formula: 


^-{aiixj+a22xl+2ai2XlX2)+blXl+b2X2^^^^^^ _ TT ^°-12blb2+ailb2 


V (ietj4 


where: 


On ~ 2 2)?— 2gXcj) + + 2gX4>), 

022 = “ + A^ + 2(j)^ + ‘2^9^4^ T — 2gX(f)), 


ai 2 = 0.21 = (A + 1 - A - (() - AgX4>)(j), 

det 24 = (^ + 2 A 2(()2 - 25 A(/> + 02((/,2 + 2gX(t))){^ + A^ + 20 ^ + 25 A 0 + (^ 2 (a 2<^2 _ 2gXcj)))- 


- {X^r + l-X^-r- 4gX(j)yr) 

Then we obtain: Coooo = 

To calculate the following integrals we differentiate Eq. 33 by the corresponding Uij : 


(33) 

(34) 

(35) 

(36) 

(37) 


Coon — 
Coiio = 


{ai2{(y - 1 ) - 4>a22 + <paii) 

(77det^)3/2VA 

{.4>oii — ai 2 ) 


Coioi — 


{4>ai2 + an) 


C3(det T)3/2((/,2 + A2 + 2£/A(/>) 1/4 ’ 


Ciooi = — 


c:3(det 24)3/2(A2<^2 + I _ 2gX(t)Yl^ 

^(y + X^ + 2gX(t)^ 


K{det ^)3/2 


rr \/A2(()2 + 1 - 2c/A(/> OuVX 

ClOlO = -- [022 - 9012 ), Ciioo = - 


3aT 


Clin = 


C(det .4)3/2 ^ K(det .4)3/2 

As the corresponding functions in the following integrals are odd: 

Coooi = Cooio = Coioo = Ciooo = Com = Cion = Cmo = Cnoi = 0 


(A:detA)3/2 


(38) 


-{(P 022 + 012 ) (39) 


(40) 


(41) 
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